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We have analyzed Coulomb drag between currents of interacting electrons in two parallel one- 
dimensional conductors of finite length L attached to external reservoirs. For strong coupling, the 
relative fluctuations of electron density in the conductors acquire energy gap M. At energies larger 
than r = const x V- exp(— LM/v-)/L + T+, where T+ is the impurity scattering rate, and for 
L > V-/M, where V- is the fluctuation velocity, the gap leads to an "ideal" drag with almost equal 
currents in the conductors. At low energies the drag is suppressed by coherent instanton tunneling, 
and the zero-temperature transconductance vanishes, indicating the Fermi liquid behavior. 
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Current drag between two parallel conductors provides 
important information about electron correlations both 
between the conductors and inside them. The domi- 
nant drag mechanism is the momentum transfer between 
electrons in the two conductors by Coulomb scattering, 
and in the case of Fermi liquids, the usual phase-space 
arguments applied to this scattering show that the lin- 
ear drag trans-conductance (and trans-resistance) should 
vanish as T 2 when temperature T goes to zero. Vanishing 
zero-temperature transconductance can be viewed then 
as a manifestation of the Fermi- liquid behavior. In this 
work, we study the low-energy behavior of the current 
drag between one-dimensional conductors of interacting 
electrons (j}-|J, where the near-ideal equality of currents 
in the two conductors was predicted || in the strong- 
coupling regime. It is shown below, however, that despite 
this "ideal" drag at large energies, the reservoirs and im- 
purity scattering suppress the drag at low energies mak- 
ing the linear zero-temperature trans- conductance zero. 
Experimentally, this behavior of the current drag can 
be studied in coupled and individually contacted one- 
dimensional AlGaAs/GaAs heterostructures [0, where 
the drag has been found recently [|||] . 

The set-up we consider consists of two identical par- 
allel ID conductors of finite length L attached to exter- 
nal reservoirs. Our description of the current drag is 
based on the sin-Gordon model. In the language of this 
model, when the conductors are infinitely long, strong 
coupling between the currents in the two conductors is 
equivalent to opening of the energy gap M in the spec- 
trum of the relative electron density fluctuations. This 
requires an interacting repulsive Tomonaga-Luttinger liq- 
uid (TLL) in the individual conductors. For conduc- 
tors of finite length L, the strong coupling is expected 
if M > Tl = V-/L, where w_ is velocity of the relative 
electron density fluctuations. Here and below e = % = 1. 
Deviations between the currents in the two conductors 
in this regime are due to tunneling of instantons of en- 
ergy M. To describe this tunneling we derive below a 
low-energy model using the Schmid's duality transfor- 



mation M, and solve it by fermionization. The main 
prediction that follows from the solution is the existence 
of a new low-energy regime of coherent multi-instanton 
tunneling that leads to a complete suppression of the 
Coulomb drag at small temperatures and voltage differ- 
ences V across the two conductors. Inclusion of a weak 
one-electron impurity scattering with small rate T + sup- 
presses the trans-conductance further, increasing the en- 
ergy range T of the coherent multi-instanton tunneling 
as T = const x exp(—LM/v-)v-/L + T + . 

Transport through a one-channel wire confined be- 
tween two reservoirs of spinless electrons and close to 
a screening gate can be modeled HQ by a ID sys- 
tem of electrons whose interaction is local and switched 
off outside the finite length of the wire. Haldane's 
bosonization pj applied to each of the two wires re- 
lates their ID electron densities pb (b — 1,2) to the 
appropriate bosonic fields <f>b as pt,(x,t) — (2kbF + 
d x 4>b{x, t)) exp{in(kbF% + ^{x, 0/2)}/27r. Here sum- 
mation runs over even n, and kbF, b = 1, 2, are the Fermi 
momenta of the two wires. Evolution of each bosonic 
field 4>b(x,t) is described with the Lagrangian of an in- 
homogeneous TLL 




where the coordinate x was scaled by inverse veloc- 
ity. The inhomogeneous interaction constant of the TLL 
model, g(x) = 1 + (g— l)(p(x), ip(x) = 9{x)9{L — x), cor- 
responds to non-interacting electrons outside the wire, 
g(x) = 1, and takes on some interacting value g(x) = 
g < 1 inside it. The fields <j>± = {<j)\ ± 02)/v / 2 related to 
the fluctuations of the total and relative electron densities 
in the wires have the Lagrangians C± of the same form, 
and the total Lagrangian is L = ^ ± C± = J2 b =i 2 
Inter-wire interaction adds another part to the total La- 
grangian, C int = - Jq dx dyU(x- y)pi(x)p 2 (y), where 
U(x — y) is short-ranged on the scale of the conductors' 
length L. Substituting bosonic expression for the densi- 
ties into Ci n t we write the total Lagrangian as 
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where f/o,i/(27r) are constants that for weak potential 
U(x) are equal to the amplitudes of the forward and 
backward scattering, and Akp = k\p — Equa- 
tion ([!]) means that the forward scattering modifies the 
constants g± of the TLL interactions and velocities of 
the modes: g± = g/ v / l± gU /ir (i.e., g + < g < gj), 
and v± = yl± gUo/ir. The inter-wire interaction also 
leads to backscattering of amplitude E F Ui/2ir 2 in the 
"—"mode, making evolution of this mode subject to the 
sin-Gordon Lagrangian. For the zero chemical potential 
/i = V-Akp — 0, a renormalization-group analysis |Tl| l 
of the uniform sin- Gordon model at energies larger than 
Tp leads to the renormalized values of the parameters 
in (|l|) for the — mode. When g_ > 1, the backscat- 
tering amplitude goes to zero if its initial value is suf- 
ficiently small, 5_ — 1 > \U\\v—/2ir. Otherwise, the 
backscattering amplitude increases and becomes of or- 
der 1 , while g_ renormalizcs to 1 /2 as the energy cut-off 
scales down to the gap Mq opening in the — mode and 
evaluated as M ~ £ f (|J7i|v_/4tt 2 (1 - ,9-)) [2(1 ~ 9 - )rl 
for 1 — g_ ^> |[/i|v_/47r 2 . A finite /x does not change 
this behavior unless /i > Mo. Since Uq > \Ui\ for re- 
alistic interactions, if g = 1, this scaling always brings 
the "-" mode into the gapless TLL regime. Therefore, 
existence of the interacting repulsive TLL's in the indi- 
vidual wires with g < 1 is a prerequisite for the opening 
of the gap in the energy spectrum of the relative density 
excitations and associated strong Coulomb drag. Assum- 
ing the gap, we construct below an effective low-energy 
model for electron transport taking into account impurity 
scattering inside the wires which introduces the term 
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to the total Lagrangian. Here A 12 (x) are the amplitudes 
of weak one-electron backscattering in the two wires. 

Duality Transformation. An effective model for ener- 
gies lower than some cut-off D 1 specified below can be 
derived from the expression for the partition function Z 
associated with the combined Lagrangians (|l|) and (||) 
following Schmid 0. Without impurities, the — and + 
modes are decoupled. Integrating out <f>- in the reser- 
voirs, we see that the "— " part of Z describes rare tun- 
neling of the massive — mode between neighboring de- 
generate vacua characterized by the quantized values of 
\[2<$>- (t, x) + 2fix = 27T771, where m is integer. Variation 
of m by ±1 corresponds to tunneling of an instanton 
(anti-instanton) through the wire. The tunneling am- 
plitude has been found as Pe~ M / TL , with the energy 



gap M — \J Mq — n 2 = Mo sin -07, from evaluation || 
of the instanton mass Mo , and by mapping onto a free 
fermionic model |L2| . The instanton calculation [O also 
gives the prefactor P = C x VI? 7 (sin 3 n7M T L ) 1 / T up to 
a constant C of order 1. A high-energy cut-off D 1 of 
the long-time asymptotics of the instanton-instanton in- 
teraction F(t) = 1ii{-\/t 2 + 1/D' 2 } created by the reser- 
voirs varies with \i from D' ~ \J MqTl for // = 0, to 
D' ~ [M /fi)T L for ft > T L . 

Weak impurity pinning does not change the form of in- 
stantons or their dynamics as long as Ep J dx\Ai t 2(x)\ <C 
M. However, it couples the ± modes, since the impurity 
Lagrangian (Q) restricted to the mth vacuum of the — 
mode creates a scattering potential for the + mode 
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The sign (— l) m of this potential varies as vacuum is 
switched from m to mil by instanton tunneling. For 
energies smaller than T' L = v+/L, the Lagrangian (j^) 
is equivalent to that of an effective point-like scatterer, 
(-l) m (2V lmp Tl/Tr)cos((/) + {T,0)/V2) 1 in the uniform 
non-interacting (g = 1) TLL, where the amplitude Vi mp 
is proportional to A± t 2 in the lowest order [ fl2"|| : Vi mp ~ 

(Ep/T^y-s+^ljdx (T,j A 3( x ) et{ ~ iyt " x ) I- Variation of 
the amplitude sign in Eq. (|^) can be accounted for by 
an auxiliary pseudospin variable with the corresponding 
Pauli matix 0-3. The energy splitting becomes an opera- 
tor o-3(2Vi mp T' L /n) cos(0 s (t, xq)/V2) acting on the pseu- 
dospin, and every (anti-)instanton tunneling reverses the 
03 values with the Pauli matrix o\. The partition func- 
tion can then be written as 
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Here So[4>+] = f dr£ + with g + = 1 is the free TLL 
Eucleadian action, and T denotes time-ordering. All 
r-integrals run from to inverse temperature 0, and 
J2j a j — 0- To have all o\ ^-matrices time-ordered in 
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(||) , we attributed each of them to a corresponding time 
t assuming that their time evolution is trivial. Noticing 
that the interaction F coincides with the pair correlator 
of some bosonic field whose evolution is described with 
the free TLL action <Sq[0_], we rewrite (||) ascribing the 
factors exp{=Fz0_(rj,O)/\/2} to the 73 instantons/anti- 
instantons: 

ZocTrv |t J £»0 + £>6i_ e -'So[0+]-5o[e-]+s'[<Ti,3,0 + ,e-]| 
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The action in (^|) is similar to a model of point scatterer 
with internal degree of freedom in the two-component 
TLL. This "dual" model (^) is equivalent to the initial 
one ([!]) at low energies. Its imaginary time evolution is 
described in the operator formalism as Z oc Tr{e~^ n } 
with a Hamiltonian 7i related to the action in (g) in a 
standard way. The model can be solved exactly through 
fermionization. 

Fermionization. From the commutation relations and 
hermiticity, the Pauli matrices can be written as a a — 
E/3, 7 ^^Uii with Majorana fermions £1,2,3 
and antisymmetrical tensor e: e 123 = 1. Point-like nature 
of the interaction in (J5J) , and appropriate time depen- 
dence of the 0_ and <j) + correlators, allow us to fermionize 
them introducing the operators: 



Here ^>±(0) are the x — values of the fermionic fields 
attributed to fluctuations of the total and relative densi- 
ties. These fields have linear dispersion relations "inher- 
ited" from the bosonic fields, with the cut-offs T' L and D', 
respectively. Substitution of these fields into the Hamil- 
tonian of the dual model (g) makes it free-electron-like 
with the interaction reduced to tunneling between the 
tp± fermions and the Majorana fermion £ 2 = £■ 

Application of voltages Vi,2 across the first and second 
wires can be described [jTjj with a gauge transformation 
0i,2 — > 01,2 — Vifit in the real-time Lagrangian (||). It is 
equivalent to transformation <j>± — > 0i,2 — V2V±t, V± = 
(Vi ± V2)/2, in the arguments of the cos-terms in the 
Lagrangians (Q) and (^|). Since each instanton tunnel- 
ing in the — mode transfers charge A(j>-/y/2n = I and 
energy — this transformation of 0_ causes a shift 
6-/V2^ 9-/V2 + V-t of the cos-argument in the real- 
time form of the action (||) . Assuming that both voltages 
are sufficiently small, \Vi^\ < Tl < M, we neglect their 
effect on the other parameters. The real-time Lagrangian 
associated with the fcrmionized Hamiltonian is: 



where the rate of impurity scattering is T + = 
2T^V 2 mp /n and the rate of the instanton tunneling 

is F_ = 2irC 2 \Jt l Mq sin 3 we- 2M / Th . The currents 
flowing through the first and second wire can be re- 
(5) lated as Ji,2 = (J+ ± J-)/V2 to the currents of 
the ± modes J± = —dt(t>±/2it. Duality Transforma- 
tion makes J_ equal to the tunneling current J- = 
-z v /r_/2ft/>±(0, i)£(t)e- iVLt - h.c] in the fcrmionic 
model (JoJ) . Meanwhile, the tunneling current of the + 
mode j+ = i^T + /2[iJj+(0,t)£(t)e lV + t - h.c] coincides 
with the backscattered current of the bosonic Lagrangian 
(||) , whose average is simply related to the average direct 
current (TJ: (J+) = V+/(V2tt) - 

Currents. Since the Lagrangian (^) is Gaussian, both 
currents can be found with the non-equilibrium Keldysh 
technique as (J_) = (\/2r_/7r) J(VL/T, T/T) and = 
(v / 2r + / 7 r)J(V+/T,r/T), where 
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In spite of the coupling between the ± modes, each of the 
currents (J±) is independent of the voltage V T applied 
to the opposite mode. At T = 0, eq. (^) gives for the 
currents through the first and the second wire: 
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The effect of the Coulomb drag is clearly seen when a 
voltage V is applied to the first (drive) wire only, V+ — 
V- = V/2. In this case, the average drag current in the 
second wire ( J 2) = — £ arctan(V~/4r) depends only on 
the total crossover energy F. The current is suppressed 
below this energy and approaches a half of the maximum 
total current agV above it, where <to = e 2 /2irh is the con- 
ductance of a ballistic ID conductor. This shows that the 
Fermi-liquid reservoirs and one-electron impurity scatter- 
ing have similar effect of suppressing the drag current. 
However, both mechanisms affect the drive wire current 
differently, ( J x ) = ^ + r ~~ r+ arctan(V/4r). Below the 
crossover, this current reaches maximum <j$V in absence 
of disorder, while the disorder suppresses it. Above the 
crossover, the current is about a half of the maximum, 
and the conductance (J\)/V takes on the fractional value 
cto/2, which is not affected by the disorder. 

In the regime linear in V±, the transport is char- 
acterised by the "diagonal" conductance Gc and 
trans-conductance Gt- To find them we use 
that J{V/T,T/T) = Vip' (1/2 + T / (irT)) / (2ttT) = 
G{T/T)V/(2T) for small V. Here ij>'(x) is the deriva- 
tive of the di-gamma function, and the function G de- 
creases monotonically from one at zero temperature, 
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G = 1 — ^(ttT/2T) 2 , to zero at large temperatures, 
G = (nT/2T) - 14C(3)(7rr/2T) 2 , where the zeta function 
C(3) w 1.2. The linear trans-conductance Gt — 0o(l — 
G)/2 is a function of T/T; it approaches 1/2 at large tem- 
peratures T ^S> T and vanishes at T <C T. The diagonal 
conductance G c = cr [1/2 + G(T/T)(T- - T + )/T] has a 
fractional large-temperature asymptotics o"o[l/2 + (T_ — 
r + )7r/(4T)]. The sign of T_ - T + determines whether 
the conductance increases or decreases with decreasing 
temperature. At zero temperature, Gc = (TqT-/T and 
varies from for T + 3> F_ to its maximum <jq in absence 
of disorder. 
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FIG. 1. Trans-resistance Rt due to the current drag be- 
tween two single-mode conductors as a fnction of temperature 
for different rates F+ of impurity scattering in the conductors. 



In experiments, the drag trans-resistance Rt is used 
typically (see, e.g. as ^ ne quantitative measure of 

strength of the Coulomb drag. We obtain Rt by in- 
version of the conductance matrix, Rt — Gt/(G c — 
G T ), and show the results in Fig. 1. Without im- 
purity scattering (r+ = 0), Gc + Gt = co and the 
trans-resistance reduces to a one-parameter function, 
Rt{T/TJ) = (1 - G{T/T-))/[2a G]. It increases with 
temperature from zero at T = first quadratically, 
R T = 7r[i(7rT/2r_) 2 - ^(7rT/2r_) 4 ] and then, above 
r_, linearly: R T = 2T/r_ + [56C(3)/tt 2 - tt] up to an 
exponentially large value Rt oc exp[M/Tz] at the up- 
per boundary of applicability of our model T ~ Tl- 
The impurity scattering changes the trans-resistance to 
R T = R T (T/T)^-[1-^F-G(T/T)]- 1 . In particular, R T is 
always enhanced by the scattering at temperatures larger 
than r as R T = R T {T/Y-) + §|C(3)r+/r_ and remains 
small below r_ (see Fig. 1). 

The results obtained in this work are valid for the fol- 
lowing hierarchy of energies, V,T <C Ij, <C M, and for 
sufficienly small amplitude of impurity scattering, Vi mp < 
M/T L (T L /E F ) a +/ 2 . They show that without impurity 
scattering, the linear trans-conductance Gt is exponen- 
tially close to (Jo/2 at T <C T < Tl. Gt should remain 



exponentially close to Co/2 even at larger temperature, 
T >Tl, but with deviations from this value now due to 
the thermal activation of instantons, Gt—ctq/2 cx eT M l T . 
With further increase of temperature beyond M, Gt is 
expected to deviate from <7o/2, decreasing as 1/T, as 
can be seen from a solution |l2] ] of a similar problem. 
When Tl approaches M (e.g., with decreasing conductor 
length) the width of the peak in the Gt(T) dependence 
around Tl reduces and this dependence smoothly goes 
over into the perturbative regime. In this case, similarly 
to the strong coupling regime, Gt oc T 2 at low temper- 
atures, T < T L , but G T oc T 4 9-" 3 at T > T L . The 
trans-conductance can even grow with T, in particular, 
Gt oc T (l^Jl^,|^] for ID Fermi liquid conductors above 
Tl- Note, however, that the Fermi liquid description of 
the drag problem is inconsistent, since it takes into ac- 
count only the intrawire part of the Coulomb interaction. 

To evaluate the minimum length — v/M$ necessary 
for observation of the strong Coulomb drag, we notice 
that in the optimum situation for observation of drag, 
the wire parameters: the width w, their separation d, 
and the distance D from each wire to its screening gate, 
are related by the set of inequalities: w < D < d ~ l/kp. 
Standard expressions for the electrostatics of inter-wire 
and intra-wire interactions show that these inequalities 
provide sufficiently small TLL interaction parameter g, 
and not-too-small backscattering amplitude Ui of the in- 
terwire interaction, while keeping the upward renormal- 
ization of <7_ caused by the forward scattering Uq to a 
minimum. We evaluate L* ~ 3/im JTij ] for w = lOnm, 
D = 20nm, d = 40nm, kp 1 = 30nm. These parameters 
are realistic for the present-day GaAs heterostructures 
and close to satisfying the inequalities. 

In conclusion, we have shown that strong Coulomb 
drag occurs between currents of repulsive TLL in two ID 
conductors at temperatures and/or voltage differences 
above some crossover energy V and below the energy gap 
M of the relative density fluctuations. The drag is sup- 
pressed below T and grows to the ideal one at T, V — > Tl- 

* On leave of absence from A.F.Ioffe Physical Technical 
Institute, 194021, St. Petersburg, Russia. 
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